Voltage-Current curves for small Josephson junction arrays 
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We compute the current voltage characteristic of a chain of identical Josephson circuits charac- 
terized by a large ratio of Josephson to charging energy that are envisioned as the implementation 
of topologically protected qubits. We show that in the limit of small coupling to the environment 
it exhibits a non-monotonous behavior with a maximum voltage followed by a parametrically large 
region where V oc 1/1. We argue that its experimental measurement provides a direct probe of the 
amplitude of the quantum transitions in constituting Josephson circuits and thus allows their full 
characterization. 



I. INTRODUCTION 

In the past years, the dramatic experimental progress 
in the design and fabrication of quantum two level sys- 
tems in various superconducting circuits^ has raised a 
hope that such solid state devices could eventually serve 
as basic logical units in a quantum computer (qubits). 
However, a very serious obstacle on this path is the ubiq- 
uitous decoherence, which in practice limits the typical 
life-time of quantum superpositions of two distinct log- 
ical states of a qubit to microseconds. This is far from 
being sufficient to satisfy the requirements for implement- 
ing quantum algorithms and providing systematic error 
correction^ 

This has motivated us to propose some alternative 
ways to design Solid-State qubits, that would be much 
less sensitive to decoherence than those presently avail- 
able. These protected qubits are finite size Josephson 
junction arrays in which interactions induce a degener- 
ate ground-state space characterized by the remarkable 
property that all the local operators induced by couplings 
to the environment act in the same way as the identity 
operator. These models fall in two classes. The first 
class is directly inspired by Kitaev's program of topo- 
logical quantum computation^ and amounts to simulat- 
ing lattice gauge theories with small finite gauge groups 
by a large Josephson junction lattice4^£ The second 
class is composed of smaller arrays with sufficiently large 
and non-Abelian symmetry groups allowing for a persis- 
tent ground-state degeneracy even in the presence of a 
noisy environment 2& All these systems share the prop- 
erty that in the classical limit for the local superconduct- 
ing phase variables (i.e. when the Josephson coupling is 
much larger than the charging energy) , the ground-state 
is highly degenerate. The residual quantum processes 
within this low energy subspace lift the classical degen- 
eracy in favor of macroscopic coherent superpositions of 
classical ground-states. The simplest example of such 
system is based on chains of rhombi (Fig. [T]) frustrated 
by magnetic field flux $ = $0/2 that ensures that in the 
classical limit each rhombus has two degenerate states^ 

Practically, it is important to be able to test these ar- 



rays and optimize their parameters in relatively simple 
experiments. In particular one needs the means to verify 
the degeneracy of the classical ground states, the pres- 
ence of the quantum coherent processes between them 
and measure their amplitude. Another important pa- 
rameter is the effective superconducting stiffness of the 
fluctuating rhombi chain. The classical degeneracy and 
chain stiffness can be probed by the experiments dis- 
cussed m£; they are currently being performed^. The 
idea is that a chain of rhombi threaded individually by 
half a superconducting flux quantum, the non-dissipative 
current is carried by charge 4e objects ; 11 ' 12 so that the 
basic flux quantum for a large closed chain of rhombi 
becomes h/(4e) instead of h/(2e) which can be directly 
observed by measuring the critical current of the loop 
made from such chain and a large Josephson junction. 

The main goal of the present paper is to discuss a prac- 
tical way to probe directly the quantum coherence associ- 
ated with these tunneling processes between macroscop- 
ically distinct classical ground-states. In principle, it is 
relatively simple to implement, since it amounts to mea- 
suring the average dc voltage generated across a finite 
Josephson junction array in the presence of a small cur- 
rent bias (i. e. this bias current has to be smaller than 
the critical current of the global system). The physi- 
cal mechanism leading to this small dissipation is very 
interesting by itself; it was orinally discussed in a sem- 
inal paper by Likharev and Zorini£ in the context of a 
single Josephson junction. Consider one element (single 
junction or a rhombus) of the chain, and denote by <f> 
the phase difference across this element. When it is dis- 
connected from the outside world, its wave-function \E' 
is 27r£-periodic in <p where £ = I for a single junction 
and £ = 1/2 for a rhombus. This reflects the quanti- 
zation of the charge on the island between the elements 
which can change by integer multiples of 2e/£. If <fi is 
totally classical, the element's energy is not sensitive to 
the choice of a quasi-periodic boundary condition of the 
form ^(c/) + 2tt() = cxp^irCq)^ (<j>) , where q represents 
the charge difference induced across the rhombus. In the 
presence of coherent quantum tunneling processes for </>, 
the energy of the element e(q) will acquire g-dependence, 
with a bandwidth directly related to the basic tunnel- 
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ing amplitude. Whereas q is constrained to be integer 
for an isolated system, it is promoted to a genuine con- 
tinuous degree of freedom when the array is coupled to 
leads and therefore to a macroscopic dissipative environ- 
ment. So, as emphasized by Likharev and Zorini&, the 
situation becomes perfectly analogous to the Bloch the- 
ory of a quantum particle in a one-dimensional periodic 
potential, where the phase <f> plays the role of the posi- 
tion, and q of the Bloch momentum. A finite bias cur- 
rent tilts the periodic potential for the phase variable, so 
that in the absence of dissipation, the dynamics of the 
phase exhibits Bloch oscillations, very similar to those 
which have been predicted^ and observe d 18 ! 19 for elec- 
trons in semi-conductor super-lattices. If the driving cur- 
rent is not too large, it is legitimate to neglect inter-band 
transitions induced by the driving field, and one obtains 
the usual spectrum of equally spaced localized levels of- 
ten called a Wannier-Stark ladder. In the presence of 
dissipation, these Wannier-Stark levels acquire a finite 
life-time, and therefore the time-evolution of the phase 
variable is characterized by a slow and uniform drift su- 
perimposed on the faster Bloch oscillations. This drift is 
translated into a finite dc voltage by the Josephson re- 
lation 2eV = H(d<f>/di). This voltage decreases with cur- 
rent until one reaches the current bias high enough to 
induce the interband transition. At this point the phase 
starts to slide down fast and the junction switches into 
a normal state. In the context of Josephson junctions 
these effects were first observed in the experiments on 
Josephson contacts with large charging energ y 20 ! 21 i 22 i 23 
and more recentl y 24 ! 25 in the semiclassical (phase) regime 
of interest to us here. Bloch oscillations in the quantron- 
ium circuit driven by a time-dependent gate voltage have 
also been recently observed.— 

This picture holds as long as the dissipation affecting 
the phase dynamics is not too strong, so that the radia- 
tive width of the Wannier-Stark levels is smaller than the 
nearest-level spacing (corresponding to phase translation 
by 27r£) that is proportional to the bias current. This 
provides a lower bound for the bias current which has 
to be compatible with the upper bound coming from the 
condition of no inter-band transitions. As we shall see, 
this requires a large real part of the external impedance 
Z u ^> Rq as seen by the element at the frequency of the 
Bloch oscillation, where the quantum resistance scale is 
Rq = h/ (4e 2 ). This condition is the most stringent in or- 
der to access experimentally the phenomenon described 
here. Note that this physical requirement is not lim- 
ited to this particular experimental situation, because 
any circuit exploiting the quantum coherence of phase 
variables, for instance for quantum information process- 
ing, has to be imbedded in an environment with a very 
large impedance in order to limit the additional quantum 
fluctuations of the phase induced by the bath. The intrin- 
sic dissipation of Josephson elements will of course add to 
the dissipation produced by external circuitry, but we ex- 
pect that in the quantum regime (i.e. with sizable phase 
fluctuations) considered here, this additional impedance 



will be of order of Rq at the superconducting transition 
temperature, and will grow exponentially below. Thus, 
the success of the proposed measurements is also a test of 
the quality of the environment for the circuits intended 
to serve as protected qubits. 

In many physical realizations has a significant fre- 
quency dependence and the condition Z^ 3> Rq is sat- 
isfied only in a finite frequency range uj max > ui > Lu m - m . 
This situation is realized, for example, when the Joseph- 
son element is decoupled from the environment by a 
long chain of larger Josephson junctions (Section fV| . 
In this case the superconducting phase fluctuations are 
suppressed at low frequencies implying that a phase co- 
herence and thus Josephson current reappears at these 
scales. The magnitude of the critical current is however 
strongly suppressed by the fluctuations at high frequen- 
cies. This behavior is reminiscent of the reappearance 
of the Josephson coupling induced by the dissipative en- 
vironment observed ir*2i. At higher energy scales fluc- 
tuations become relevant, the phase exhibits Bloch os- 
cillations resulting in the insulating behavior described 
above. Thus, in this setup one expects a large hierarchy 
of scales: at very low currents one observes a very small 
Josephson current, at larger currents an almost insulat- 
ing behavior and finally a switching into the normal state 
at largest currents. 

In the case of a chain of identical elements, the total 
dc voltage is additive, but Bloch oscillation of different 
elements might happen either in phase or in antiphase. 
In the former case the ac voltages add increasing the 
dissipation in the external circuitry; while in the latter 
case the dissipation is low and the individual elements 
get more decoupled from the environment. As we show in 
Section IIIII a small intrinsic dissipation of the individual 
elements is crucial to ensure the antiphase scenario. 

This paper is organized as follows. In section [Til we 
present a semi-classical treatment of the voltage versus 
bias current curves for a single Josephson element. We 
show that this gives an accurate way to measure the effec- 
tive dispersion relation e{q) of this element, which fully 
characterizes its quantum transition amplitude. Further, 
we show that application of the ac voltage provides a 
direct probe of the periodicity {2tt versus 7r) of each ele- 
ment. In Section IIIII we consider the chain of these ele- 
ments and show that under realistic assumptions about 
the dynamics of individual elements, it provides much 
more efficient decoupling from the environment. Sec- 
tion IIVI focusses on the dispersion relation expected in 
a practically important case of a fully frustrated rhom- 
bus which is the building block for the protected arrays 
considered before^ In this case, the band structure has 
been determined by numerical diagonalizations of the 
quantum Hamiltonian. An important result of this anal- 
ysis is that even in the presence of relatively large quan- 
tum fluctuations, the effective band structure is always 
well approximated by a simple cosine expression. Finally, 
in section [V] we discuss the conditions for the experimen- 
tal implementation of this measurement procedure and 
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the full V(I) characteristics expected in realistic setup. 
After a Conclusion section, an Appendix presents a full 
quantum mechanical derivation of the dc voltage when 
the bias current is small enough so that inter-band tran- 
sitions can be neglected, and large enough so that the 
level decay rate can simply be estimated from Fermi's 
golden rule. 



II. SEMI-CLASSICAL EQUATIONS FOR A 
SINGLE JOSEPHSON ELEMENT 



Let us consider the system depicted on Fig. ([I]). In 
the absence of the current source, the energy of the one 
dimensional chain of N Josephson elements is a 2ir( pe- 
riodic function of the phase difference <fi across the chain. 
The current source is destroying this periodicity by in- 
troducing the additional term —h(I/2e)(f) in the system's 
Hamiltonian. Because </> is equal to the sum of phase 
differences across all the individual elements, it seems 
that the voltage generated by the chain is iV times the 
voltage of a chain reduced to a single element. This is, 
however, not the case: the individual elements are cou- 
pled by the common load, and furthermore, as we show 
in the next section, their collective behavior is sensitive 
to the details of the single element dynamics. In this sec- 
tion, we consider the case of a single Josephson element 
(N = 1), rederive the results of Likharev and Zorini& for 
single Josephson contact and generalize them for more 
complicated structures such as rhombus and give ana- 
lytic equations convenient for data comparison. 

The dynamics of a single Josephson contact is analo- 
gous to the motion of a quantum particle (with a charge 
e) in a one-dimensional periodic potential (with period 
a) in the presence of a static and uniform force F, the 
phase-difference <f> playing the role of the spatial coordi- 
nate x of the particle^ In the limit of a weak external 
force, it is natural to start by computing the band struc- 
ture e n (k) for k in the first Brillouin zone [— 7r/a, 7r/a], 
n being the band label. A first natural approximation is 
to neglect interband transitions induced by the driving 
field. This is possible provided the Wannier-Stark energy 
gap Ab = Fa is smaller than the typical band gap A in 
zero external field. As long as Ab is also smaller than 
the typical bandwidth W, the stationary states of the 
Schrodinger equation spread over many (roughly PF/Ab) 
periods, so we may ignore the discretization (i.e. one 
quantum state per energy band per spacial period) im- 
posed by the projection onto a given band. We may 
therefore construct wave-packets whose spacial extension 
Ax satisfies a <C Ax <§; aWj Ab , and the center of such a 
wave-packet evolves according to the semi-classical equa- 
tions: 



dx 
It 
dk 
It 



1 de n (k) 
h dk 
1 



h 



-F 



(1) 
(2) 



In the presence of dissipation, the second equation is 
modified according to: 

where m* is the effective mass of the particle in the n-th 
band and r is the momentum relaxation time introduced 
by the dissipation. 
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FIG. 1: The experimental setup discussed in this paper: a 
chain of identical building blocks represented by shaded rect- 
angle that are biased by the external current source charac- 
terized by the impedance Z{lu). The internal structure of the 
block that is considered in more detail in the following sec- 
tions is either a rhombus (4 junction loop) frustrated by half 
flux quantum, or a single Josephson junction but the the re- 
sults of the section [TT] can be applied to any circuit of this 
form provided that the junctions in the elementary building 
blocks are in the phase regime, i.e. E,j 2> Ec- 

In the context of a Josephson circuit, we have to diago- 
nalize the Hamiltonian describing the array as a function 
of the pseudo-charge q associated with the 2tt( periodic 
phase variable <j). The quantity q controls the periodic 
boundary condition imposed on (f>, namely the system's 
wave-function is multiplied by exp(i2irq) when (f> is in- 
creased by 2tt(. From this phase-factor, we see that 
the corresponding Brillouin zone for q is the interval 
[—1/2, 1/2]. For a simple Josephson contact (£ = 1), the 
fixed value of q means that the total number of Cooper 
pairs on the site carrying the phase <j) is equal to q plus an 
arbitrary integer. For a doubly periodic element, such as 
rhombus (£ = 1/2), charge is counted in the units of 4e. 
To simplify the notations we assume usual 27r periodicity 
(C = 1) in this and the following Sections and restore 
the £- factors in Sections IIV1 fVl From the band structure 
e n(q), we may write the semi-classical equations of mo- 
tion in the presence of the bias current / and the outer 
impedance Z as: 

# 1 de n (q) 

dt h dq 

dq I Zq d(f> 

~dt ~ 2e ~~Z~~dt 
where we used the Josephson relation for the voltage drop 
V across the Josephson element as V = (h/ 2e) (dtp / dt) 
and defined Z Q = h/{Ae 2 ). 



(4) 
(5) 



4 



This semi-classical model exhibits two different 
regimes. Let us denote by w max the maximum value of 
the "group velocity" \de n (q)/(hdq)\. If the driving cur- 
rent is small (7 < I c = 2eu) max ZQ/Z), it is easy to see 
that after a short transient, the system reaches a station- 
ary state where q is constant and: 



dt 



I Z 
Yz^ 



(G) 



that is: V = ZI. Thus, at / < J c the current flows en- 
tirely through the external impedance, i.e. the Joseph- 
son elements become effectively insulating due to quan- 
tum phase fluctuations. Indeed, a Bloch state writ- 
ten in the phase reprentation corresponds to a fixed 
value of the pseudo-charge q and non-zero dc voltage 
(l/2e)(de n /dq). Note that the measurement of the maxi- 
mal value V c of the voltage on this linear branch directly 
probes the spectrum of an individual Josephson block, 
because V c = Hui max /2e 

At stronger driving (I > I c ), it is no longer possible 
to find a stationary solution for q. The system enters 
therefore a regime of Bloch oscillations. In the absence of 
dissipation (Z/Zq — > oo), the motion is periodic in time 
for both <p and q. A small but finite dissipation preserves 
the periodicity in q, but induces an average drift in <fi or 
equivalently a finite dc voltage. To see this, we first note 
that the above equations of motion imply: 



dq 
~dt 



I 

Ye 



Zq 1 de n 
Z h dq 



(7) 



Since the right-hand side is a periodic function of q with 
period 1, q(t) is periodic with the period T(I) given by: 



,1/2 

T(I)= / f(q)dq 

J -1/2 



with 



/(?) 



J_ _ Zq 1 de n 
2e Z h dq 



(8) 



(9) 



On the other hand, the instantaneous dissipated power 
reads: 



d ( , , hi 

dtV n{q) -Ye' 



Z V dt ' 



(10) 



Because q(t) is periodic, averaging this expression over 
one period gives: 



t dtj>. _ 2eZ Q d(j> , 2 . 
"dt'~ I Z {{ dt> ' 



or, equivalently: 



(ii) 



(12) 



Using the equations of motion, we get more explicitely: 



(V) 



J\%f{q)dq 



(13) 



Here we emphasized by the subscript that Z^ might have 
some frequency dependence. As we show in Appendix, 
the dissipation actually occurs at the frequency of Bloch 
oscillations that becomes u>b = 2nl/2e in the limit of 
large currents. In the limit of large currents, I 3> I c , 
(that can be achieved for large impedances) we may ap- 
proximate f(q) by a constant, so the voltage is given by 
the simpler expression: 



(V(I » 7 C )) = 



4e 2 Z 



1 f 1 ^ 2 de 
Zu>I J -i/2 dq ^ d<1 



(14) 



On the other hand, when / approaches I c from above, 
Bloch oscillations become very slow and f{q) is strongly 
peaked in the vicinity of the maximum of the group ve- 
locity. Since this velocity is in general a smooth function 
of q, we get in this limit for the maximal dc voltage: 



V r = 



; 2 
max 



4e 2 Z / c 



— ZqI c 



(15) 




FIG. 2: Typical I — V curve of a single Josephson element 
measured by a circuits shown in Fig. 1 

In the simplest case of a purely harmonic dispersion, 
e(q) = 2w cos2Trq, the maximal voltage V c — Airw/(2e). 
If one can further neglect the frequency dependency of 
Z, the V(J) can be computed analytically: 



(V) = ZI 
(V) = ZI C 



KIc 

Ic 



i+Vp^i! 



I>Ir 



(16) 
(17) 



We show this dependence in Fig. [2] This expression (|16ll . 
(fT"7f is related to the known result for Z <C Zq 1 -^ by 
the duality^ transformation: 



V 
I 



I: 

v, 

1 

z' 



The semi-classical approximation is valid when the os- 
cillation amplitude of the superconducting phase is much 
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larger than 2ir, which allows the formation of the semi- 
classical wave-packets. When / is much larger than I c , 
this oscillation amplitude is equal to 2eW/hI, where W is 
the total band- width of e n (q) . This condition also ensures 
that the work done by the current source when the phase 
increases by 2tt is much smaller than the band-width. In 
order to observe the region of negative differential resis- 
tance, corresponding to the regime of Bloch oscillations, 
we require therefore that: 



2irhl c 
2e 



< W 



2eK 



(18) 



where the last equality becomes exact in the case of a 
purely harmonic dispersion. This translates into: 



Z> R Q . 



(19) 



For large currents one can compute dc voltage di- 
rectly by using the golden rule (without semiclassics) ; 
we present the results in Appendix [X] The result is con- 
sistent with the large I limit of Eq. ([17]), (V) = V C 2 /2ZI. 
Deep in the classical regime (Ej ^> Ec), the bandwidth 
and the generated voltage become exponentially small. 
In this regime the bandwidth is much smaller than the 
energy gaps, so these formulas are applicable (asuming 
(TT9")) is satisfied) until the splitting between Wannier- 
Stark levels becomes equal to the first energy gap given 
by the Josephson plasma frequency, i.e. for / < eujj/n. 
Upon a further increase of the driving current in this 
regime the generating voltage experiences resonant in- 
crease for each splitting that is equal to the energy gap: 
Ik = e(Ek—Eo) I tt. Physically, at these currents the phase 
slips are rare events that lead to the excitation of the 
higher levels at a new phase value that are followed by 
their fast relaxation. At very large energies, the band- 
width of these levels becomes larger than their decay rate 
due to relaxation, (Rq/Z)Ec- At these driving currents, 
the system starts to generate large voltage and switches 
to a normal state. At a very large Ej this happens at 
the driving currents very close to the Josephson criti- 
cal current 2eEj, but in a numerically wide regime of 
100 > Ej /Ec > 10 the generated voltage at low curents 
is exponentially small but switching to the normal state 
occurs at significantly smaller currents than 2eEj. 

In the intermediate regime where Ej and Ec are com- 
parable, we expect a band-width comparable to energy 
gaps so that the range of application of the quantum 
derivation is not much larger than the one for the semi- 
classical approach. 

Negative differential resistance associated to Bloch os- 
cillations has been predicted long ago^ and observed 
experimentally 2 ^ in the context of semi-conductor su- 
perlattices. For Josephson junctions in the cross-over 
regime (Ej/Ec — 1), a negative differential resistance 
has been observed in a very high impedance environ- 
ment^ in good agreement with earlier theoretical pre- 
dictions. 30 More recently, the I — V curve of the type 
shown on Fig [2] have been reported on a junction with 
a ratio Ej/Ec = 4.5 25 . These experiments show good 



agreement with a calculation which takes into account 
the noise due to residual thermal fluctuations in the 
resistor i^L 

Although the above results allows the extraction of the 
band structure of an individual Josephson block from the 
measurement of dc I — V curves, the interpretation of ac- 
tual data may be complicated by frequency dependence 
of the external impedance Z u . Additional information 
independent on Z u can be obtained from measuring the 
dc V(I) characteristics in the circuit driven by an ad- 
ditional ac current. In this situation, the semi-classical 
equations of motion become: 



dt 

dq 

dt 



1 de n (q) 
h dq 

I + V cos(wi) 
2~e 



Zq d<f) 
~Y~dl 



(20) 
(21) 



A small ac driving amplitude /' strongly affects the V(I) 
curve only in the vicinity of resonances where tiuj^Ir) = 
muj, with m and n integers. The largest deviation occurs 
for m = n = 1. Furthermore, for I' <C I the terms 
with m > 1 are parametrically small in V jl while for 
/ 3> I c the terms with n > 1 are parametrically small 
in I c / 1. Experimental determination of the resonance 
current, Ir, would allow a direct measurement of the 
Bloch oscillation frequency and thus the periodicity of the 
phase potential (see next Section). Observation of these 
mode locking properties have in fact provided the first 
experimental evidence of Bloch oscillations in a single 
Josephson junction i 20 i 21 

We now calculate the shape of V(I) curve in the vicin- 
ity of m = n — 1 point when both I' <C I and I ^> I c . We 
denote by 4> (t) and qo(t) the time-dependent solutions 
of the equations ai I — Ir, in the absence of ac driving 
current. We shall look for solutions which remain close 
to <^ (i) and qo(t) at all times and expand them in small 
deviations <p 1 = 4> ~ 4>o> Qi = 1 ~ <Zo- We can always 
assume that q\ has no Fourier component at zero fre- 
quency because such component can be eliminated by a 
time translation applied to qo. The equations for 4>i,qi 
become 



d(p 1 1 

~dt = r Mqi 

dqi 



(22) 



_ I-Ir + I' cos(t^) Z Q d<p 1 

dt 2e Z dt { ' 

Because the main component of d ^2"^ oscillates with 
frequency u> and q\ has no dc component, the average 
value of the voltage -J^- is due to the part of q\ that os- 
cillates with the same frequency, qi^ — I'/(2eui) sin(wt). 
Because qo = LL>(t — t )+x{u(t — t )) where \(t) is a small 
periodic function, the first equation implies that 

< Itt > = ^ e "( w ^- t o)) sin M)) 



1 



cos(2irq)dq 
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The deviation qi remains small only if the constant parts 
cancel each other in the right hand side of the equation 
(|23l). This implies 



#1 

<1F >= 



z I-I R < 1 
Zpi 2e h 



,(g) cos(2irq)dq (24) 



We conclude that in the near vicinity of the reso- 
nances the increase of the current does not lead to addi- 
tional current through the Joscphson circuit, so the re- 
lation between current and voltage becomes linear again 
SV = Z5I. In other words, the Josephson circuit be- 
comes insulating with respect to current increments. The 
width of this region (in voltage) is directly related to the 
first moment of the energy spectrum of the Josephson 
block providing one with the direct experimental probe 
of this quantity. In particular, a Josephson element such 
as rhombus in a magnetic flux somewhat different from 
$o/2 displays a phase periodicity 2ir but a very strong 
deviations from a simple cos 2nq spectrum that will man- 
ifest themselves in first moment of the spectrum. Note 
finally, that the discussion above assumes that the ex- 
ternal impedance Z u has no resonances in the important 
frequency range. The presence of such resonances will 
modify significantly the observed V(I) curves because it 
would provide an efficient mechanism for the dissipation 
of Bloch (or Josephson) oscillations at this frequency. 



III. CHAIN OF JOSEPHSON ELEMENTS 

We shall first consider the simplest example of a two- 
element chain, because it captures the essential physics. 
This chain is characterized by two phase differences ((fri 
and 4> 2 ) and two pseudo-charges (q± and g 2 ). The equa- 
tions of motion for the pseudo-charges ([5]) implies that 
the charge difference q\ — qi is constant, because the cur- 
rents flowing through these elements are equal, and thus 
the right-hand sides of the evolution equations ([5]) are 
identical. Because of this conservation law, even the long- 
term physical properties depend on the initial conditions. 
Similar problems have already been discussed in the con- 
text of a chain of Joscphson junctions driven by a current 
larger than the critical currenti 32 ' 33 i 34 ' 35 This unphysical 
behavior disappears if we take into account the dissipa- 
tion associated with individual elements. Physically, it 
might be due to stray charges, two-level systems, quasi- 
particles, phonon emission, etc . 36 ' 37 A convenient model 
for this dissipation is to consider an additional resistor 
in parallel with each junction. For the sake of simplicity, 
we assume that each element has a low energy band with 
a simple cosine form. This physics is summarized by the 
equations: 



Qi 



j = Attw sin 2itqj 

2 e \ 9(.7^ n Op U , 



2eZ 



2eRj 



(25) 
(26) 



Eliminating the phases gives: 

(q\ + Qi sin 2irqi) = v ^~( sm + sm 2 7ri ?2) 

(<?2 + ^2 sin27rg 2 ) = v — ^ (sin27rgi + sin27rg 2 ) 



where 
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Ye 
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Here we allowed for different effective resistances asso- 
ciated with each element because this has an important 
effect on their dynamics. Indeed the difference between 
the currents flowing through the resistors changes the 
charge accumulated at the middle island and therefore 
violates the conservation law mentioned before. Using 
the notations <5f2 = (f2 2 — fii)/2 and q± — (g 2 ± <7i)/2, 
we have: 



gl + f2 sin 27r<2_ cos 2?tq + + 
+SQ cos 27rq_ sin 2?tq + 
+ (i/q + f2) sin 27r<2 + cos 2irq_ — 
—SD, sin 2 7177 _ cos 2itq + 



= 



(27) 
(28) 



Significant quantum fluctuations imply that internal re- 
sistance of the element R ~ Zq for individual elements 
at T < Tc; at lower temperature it grows exponen- 
tially. Thus, in a realistic case i? 3> Z which implies that 
f?i -C v. In the insulating regime the equations (|27H28|) 
have stable stationary solution (y§ + f2) sin 2nq + = v. 
q_ = 0. This solution exists for (^o + 0) < v , i.e. if 
the voltage drop across both junctions does not exceed 
V c = 8irw j (2e) . The conducting regime occurs when 
v > (2/0 + 0); to simplify the analytic calculations we 
assume that v 3> vq- This allows to solve the equations 
(I27H28[) by iterations in all non-linear terms. In the ab- 



sence of non-linearity q^ 



vt 



const: the first 



iteration gives periodic corrections oc cos 2irvt. Averag- 
ing the result of the second order iteration over the period 
we get 



lq'-) = \v cos 2 2tt(j_ + 20] 

2v J 



(29) 



The second term in the right hand side of this equation 
is much smaller than the first if fl <C vq. In its absence 
the dynamics of g_ has fixed points at cos27rq_ = 0. At 
these fixed points the periodic potentials generated by in- 
dividual elements cancel each other and the dissipation in 
external circuitry (which is proportional to cos 2 (27rq_)) 
is strictly zero. In a general case the equation (f2"9"| has 
solution 

1 

cos 2 (27rq_) = 



^±ptan 2 



im^2VL(y Q + 2srjt 



7 



that corresponds to the short bursts of dissipation in 
external circuitry that occur with low frequency = 
^5fly / 2fl(iy + 2Q). The average value of cos 2 (2-7r<7_) 



< cos 2 (27r<7_) >= 



1 



1 



-2S2 



2Vt 



20 



i/ + 2ft 



is small implying that the effective dissipation introduced 
by the external circuitry is strongly suppressed because 
the pseudocharge oscillations on different elements al- 
most cancel each other. The effective impedance of the 
load seen by individual junction is strongly increased: 



V0 + 2Q 



(30) 



Similar to a single element case discussed in the previous 
Section, an additional dissipation in the external circuit 
implies dc current across the Josephson chain 



V 



v c f T I » I c = v c /z eff 



We conclude that a chain of Josephson elements has a 
current-voltage characteristics similar to the one of the 
single element with one important difference: the effec- 
tive impedance of the external circuitry is strongly en- 
hanced by the antiphase locking of the individual Joseph- 
son elements. In particular, it means that the condition 
Z Rq is much easier to satisfy for the chain of the 
elements than for a single element. The analytical equa- 
tions derived here describe the chain of two elements but 
it seems likely that similar suppression of the dissipation 
should occur in longer chains. 

To substantiate this claim, lets us generalize the av- 
eraging method which led to Eq. (|29p for N = 2. The 
coupled equations of motion read: 



Qj shi27rqj = v — ^ sin27r<2/ £ (31) 



fe=l 



To second order in f2j and z/ , the averaged equations of 
motion are: 

2 N 



{q) = _^L_^°y 

w/ 2v Ais ^ 



k=l 



of — ' 



k.l 



4v 



N 



^cos(27r(q i - q k )) (32) 



fe=l 



This set of coupled equations is similar to the Kuramoto 
model for coupled rotors 3 ^ defined as: 



K 

ii = ^] ~ J^^2 sili ( 27r (l3 ~ Ik) +ot) (33) 
fc=i 



N 



The equation of motion ([33)1 exhibits synchronisation of 
a finite fraction of the rotors only for K > JC c (a) i 39 i 40 
The last term in Eq. ([32| is equivalent to the interac- 
tion term of Kuramoto model with a = Tt/2. The ad- 
ditional (third) term in the model (|32|) is the same for 
all oscillators, it is thus not correlated with individual 
qj and thus can not directly lead to their synchroni- 
sation. Remarkably, it turns out that for model (|33|) 
K c (a = it/2) = 39 i 40 , suggesting that in our case, syn- 
chronization never occurs on a macroscopic scale. Note 
that the coupling K arising from Eq. (f3"2"|) is not only j- 
dependent, but it is also proportional to N. This could 
present a problem in the infinite N limit, but should not 
present a problem in a finite system. It is striking to see 
that a = 7r/2 is the value for which synchronization is 
the most difficult. 



IV. ENERGY BANDS FOR A FULLY 
FRUSTRATED JOSEPHSON RHOMBUS 



In order to apply general results of the previous sec- 
tion to the physical chains made of Josephson junctions 
or more complicated Josephson circuits we need to com- 
pute the spectrum of these systems as a function of the 
pseudocharge q conjugated to the phase across these ele- 
ments. In all cases the superconducting phase in Joseph- 
son devices fluctuates weakly near some classical value 
O where the Josephson energy has a minimum in the 
limit Ej/Ec 3> 1. In the vicinity of the minimum, the 
phase Hamiltonian is H = —AEc-^ + \E" {4>q)(4i~ 4> Q ) 2 , 
so a higher energy state of the individual element (at a 
fixed q) can be approximated by one of the oscillator 
E n = (n + where the Josephson plasma frequency 

The Josephson en- 



uj = ^8E"(<p )E c w V&EjEc. 
ergy is periodic in the phase with the period 2tt but the 
amplitude of the transitions between these minima is ex- 
ponentially small: 

w = ahjjjiEj/Ec) 1 ^ exp(-cy/Ej/E c ) 

where a, c ~ 1. In this limit one can neglect the contribu- 
tion of the excited states (separated by a large gap lu.j) 
to the lower band, so the low energy spectrum acquires 
a simple form e(q) = 2w cos2?:q. The numerical coeffi- 
cients c, a in the formulae for the transition amplitude 
depend on the element construction. For a single junc- 
tion a s = 8 2 1 / 4 /y / 7r , c s — \/8 while for the rhombus 
a r w 4.0 , c r ~ 1.61. In case of the rhombus in mag- 
netic field with flux $0/2 the Hamiltonian is periodic in 
phase with period n provided that the rhombus is sym- 
metric along its horizontal axis: indeed in this case the 
combination of the time reversal symmetry and reflec- 
tion ensures that the Josephson energy has a minimum 



for 



±7r/2. Thus, in this case the period in q dou- 



bles and the low energy band becomes e(q) = 2wcosirq. 
The maximal voltage generated by the chain of TV such 
elements at / = I c = {&irC,ew / K)(Zq/ Z) is 



8 



V r = N 



2e 



(34) 



The voltage generated at larger currents depend on the 
collective behaviour of the elements in the chain. For a 
single element it is simply 



(V(I)) = 



(2tt<» 2 



e 2 Z u I+^P^P C ' 



For more than one element the total volage is sufficiently 
reduced due to the antiphase correlations. Generally, one 
expects that 



(V(I)) = N 



(2ttCw) 5 



1 



zeff 



(35) 



where is the effective impedance of the environment 
affecting each Josephson element which is generally much 
larger than its 'bare' impedance Z w . For two elements the 
exact solution (see previous Section) gives Zjff « V RZ 
that shows th e incr ease of the effective impedance by a 
large factor \J R/Z. We expect that a similar enhance- 
ment factor appears for all N > 2. Finallly, For I < I c , 
the system is ohmic with: 



(V(I)) = Z I 



(36) 



As discussed in Section [TTJ application of a small addi- 
tional ac voltage produces features on the current-voltage 
characteristics for the currents that produce Bloch oscil- 
lation with frequencies commensurate with the frequency 
of the applied ac field cub — 2n(^I /2e = (m/n)ui. At these 
currents the system becomes insulating with respect to 
current increments, the largest such feature appears at 
m = n = 1 that allows a direct measurement of the 
Josephson element periodicity. 

For smaller Ej/Ec ~ 1 the quasiclassical formulas 
for the transition amplitudes do not work and one has 
to perform the numerical diagonalization of the quan- 
tum system in order to find its actual spectrum. As 
Ej/Ec — > f the higher energy band approaches the 
low energy band and the dispersion of the latter de- 
viates from the simple cosine form shown in Figure [3] 
These deviations, lead to higher harmonics in the dis- 
persion: e(q) = 2wcos2n(q + 2w' cos4n(q and change 
the equations (|34l35p . Our numerical diagonalization of 
a single rhombus shows, however, that even at relatively 
small Ej/Ec ~ 1 the second harmonics w' does not ex- 
ceed 0.15w, so its additional contribution to the voltage 
current characteristic (oc w' 2 ) can always be neglected. 
Thus, in the whole range of Ej/Ec > 1 the voltage cur- 
rent characteristic is given by Eqs. (|34I35[) where the ef- 
fective value of transition amplitude t can be found from 
the band width W = E x - E = 4w plotted in Fig. H 
For comparison we show the variation of the lower band 
width for a single junction in Fig. [5] 



E 2 /E c 
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FIG. 3: Spectrum of a single rhombus biased by magnetic flux 
$ = $o/2. The upper pane shows the bands of the rhombus 
characterized by Josephson eneergy Ej/Ec = 4 as a function 
of bias charge, q. The two lower levels are fitted by the first 
two harmonics (dashed line), the coefficient tv' of the second 
harmonics is w' = O.lw. One observes period doubling of the 
first two states that reflects the symmetries of the rhombus 
frustrated by a half flux quantum. The second excited level is 
doubly degenerate that makes its period doubling difficult to 
observe. Physically, these states correspond to an excitation 
localized on the upper or lower arms of the rhombus. The 
lower pane shows the dependence of the gaps for q — as a 
function of Ej/Ec- Because higher order harmonics are very 
small for all Ej/Ec > 1, the gap E\ — Eo coincides with 4w 
where w is the tunneling amplitude between the two classical 
ground states. 



V. PHYSICAL IMPLEMENTATIONS 

Generally, the effects described in the previous sections 
can be observed if the environment does not affect much 
the quantum fluctuations of individual elements and the 
resulting quasiclassical equations of motion. These physi- 
cal requirements translate into different conditions on the 
impedance of the environment at different frequencies. 

We begin with the quantum dynamics of the individual 
elements. The effect of the leads impedance on it can 
be taken into account by adding the appropriate current 
term to the phase equation of motion before projecting 
on a low energy band and requiring that their effect on 
the phase dynamics is small at the relevant frequencies. 
For instance, for a single junction 



/ I d 2 ^ 



Zq 

Zu: 



dt 
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FIG. 4: Band width W = 4w of a single Josephson junction 



The characteristic frequency of the quantum fluctuations 
responsibe for the tunneling of a single element is Joseph- 
son plasma frequency, ujj= y/8EjE c , so the first condi- 
tion implies that 



\Z(uj)\ » y/E c /EjZQ 



(37) 



For a typical uij/2tt ~ 10GHz, the impedance of a sim- 
ple superconducting lead of the length ~ 1cm is smaller 
than Zq and the condition (|37p is not satisfied. The 
situation is changed if the Josephson elements are decou- 
pled from the leads by a large resistance or by a chain 



of M ^> 1 large junctions with y Ej/E c 1 that has 
no quantum tunneling transitions of their own (the am- 



plitude of such transitions is oc eicp(—y8Ej/E c ). As- 
suming that elements of this chain have no direct ca- 
pacitive coupling to the ground (M 2 Co <C C), the chain 

has an impedance Z = \J 8E C / EjMZq at the relevant 
frequencies, so a realistic chain with M ~ 50 junctions 



and y8Ej/E c ~ 10 provides the contribution to the 
impedance needed to satisfy (f3"7| . 

Similar decoupling from the leads of the individual el- 
ements can be achieved by a sufficiently long chain of 
similar Josephson elements, e.g. rhombi. Consider a 
long (N 3> 1) chain of similar elements connected to 
the leads characterized by a large but finite capacitance 
C g ^> C. For a short chain the tunneling of a single 
element changes the phase on the leads resulting in a 
huge action of the tunneling process. However, in a long 
chain of junctions, a tunneling of individual rhombus may 
be compensated by a simultaneous change of the phases 
8<j)/N of the remaining rhombi, and subsequent relax- 
ation of 5<j) from its initial value 7r towards the equilib- 
rium value which is zero. For N ^> 1, this later process 
can be treated within the Gaussian approximation, with 
the Lagrangian (in imaginary time) : 



L = 



1 -2 



l6E n 



2N 



--E, 



(38) 



where E g — e 2 / (2C g ). So the total action involved in the 



relaxation is: S = ^ f ^ 



1/2 



than unity this relaxation has strictly no effect on the 
tunneling amplitude of the individual rhombus. 

We now turn to the constraints imposed by the qua- 
siclassical equations of motion. The solution of these 
equations shows oscillation at the Bloch frequency that 
is uib = 2ir(I/(2e) for large currents and approaches 
zero near the I c . Thus, for a single Josephson ele- 
ment the quasiclassical equations of motion are valid if 
Re(R,Q / Z (ujb)) <C 1 ■ A realistic energy band for a 
Josephson element, W ~ 0.3K and Z/Zq ~ 100 cor- 
respond to Bloch frequency lob /27r ~ 0.1GHz . In this 
frequency range a typical lead gives a capacitive contri- 
bution to the dynamics. The condition that it does not 
affect significantly the equations of motion implies that 
the lead capacitance C < 10 fF. As discussed in Section 
pill) the individual elements in a short chain oscillate in 
antiphase decr easin g the effective coupling to the leads 
by a factor J R/Z where R is the intrinsic resistance 
of the contact. This factor can easily reach 10 2 at suf- 
ficiently low temperatures making much less restrictive 
the condition on the lead capacitance. 

Large but finite impedance of the environment 
Re(RQ / Z (ujb)) < 1 modifes the observed current- voltage 
characteristics qualitatively, specially in the limit of very 
small driving current. When / vanishes, and with infi- 
nite external impedance, the wave function of the phase 
variable is completely extended, with the form of a Bloch 
state, and the pseudo-charge q is a good quantum num- 
ber. As discussed at the end of Sec. HH the system be- 
haves as a capacitor. But when the external impedance 
is finite, charge fluctuations appear, which in the dual 
description means that quantum phase fluctuations are 
no longer unbounded. To be specific, consider a realistic 
example of N rhombi chain (or two ordinary junctions) 
attached to the leads with Z{uj) — Zq in a broad but fi- 
nite frequency interval w m i n < w < cj max and decreases as 
Z(uj) = Z (u max /uj) for lo > w max , Z(oj) = Z (uj/uj mhl ) 
for u> < w m in- Such Z(u>) is realized in a long chain 
of M Josephson junctions between islands with a finite 
capactive coupling to the ground Co: w m ax = and 
^min = {y/C/Co/M)bjj. The effective action describ- 
ing the phase dynamics across the chain has contribu- 
tions from the tunneling of individual rhombi and from 
impedance of the chain 



iujZr 



8n 2 CNw Z(uj) 



If this action S is less 



Here the first term describes the effect of the tunneling of 
the Josephson element between its quasiclassical minima 
which we approximate by a single tunneling amplitude 
w resulting in a spectrum e(q) = — 2w cos 2i:C,q that in a 
Gaussian approximation becomes e(q) = 47r 2 C 2 uig 2 . This 
approximation is justified by the fact that, as we show 
below, the main effect of the phase fluctuations comes 
from the broad frequency range where the action is dom- 
inated by the second term while the first serves only as 
a cutoff of the logarithmical divergence. Its precise form 
is therefore largely irrelevant. 
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This action leads to a large but finite phase fluctuations 



<^ 2 > = 



du> 



2lT cj 2 i iuZ Q 
8tt 2 C 2 Nw " r Z(ui) 



Z ^ min(q; max X nax ) 



where oj' max — S-k 2 ( 2 Nw(Zq/Z ). These fluctuations 
are only logarithmically large, so they result in a finite 
renormalization of the Josephson energy of the rhombi 
chain and the corresponding critical current. In the ab- 
sence of such renormalization the Josephson energy of a 
finite chain of elements can be approximated by the lead- 
ing harmonics E(<fi) = — Eq cos(</>/£) with Eq ~ Ej for 
N ~ 1 and Ej > E c . Renormalization by fluctuations 
replaces Eq by 



Ej 



1 



cxp(--<^ 2 ))i? 



minfuwXmax) 



Eo 



In the limit of w m i n — > or Z$ — > oo the phase fluc- 
tuations renormalize Josephson energy to zero. But for 
realistic parameters this suppression of Josephson energy 
is finite which thus results in a small but non-zero value of 
the critical current. In this situation the current-voltage 
characterictics sketched in Fig. Q] is modified for very 
small values of currents and voltages: instead of insulat- 
ing regime at very low currents and voltages one should 
observe a very small supercurrent (En/2e) followed by a 
small voltage step as shown in Fig. [2] by a dashed line. 
As is clear from the above discussion the value of the 
resulting critical current is controlled by the phase fluc- 
tuations at low oj <C w mal ; these frequencies are much 
smaller than the typical internal frequencies of a chain of 
Josephson elements which can be thus lumped together 
into an effective object characterized by the bare Joseph- 
son energy E(<p) and transition amplitude between its 
minima w. We thus expect the same qualitative behav- 
ior for a small chain of Josephson elements as for a single 
element at low currents. 



VI. CONCLUSION 

The main results of the present work are the expres- 
sions (|34t , (|35)) for the I-V curves of a chain of N identi- 
cal basic Josephson circuits. They are derived within the 
assumption that the Josephson coupling is much larger 
than the charging energy, but in fact, the numerical cal- 
culations show that they remain very accurate even if 
Ej ~ Ec- These equations predict a maximum dc volt- 
age when I = I c and V(I) oc l/I for I ^> I c . The anoma- 
lous V versus / dependence exhibited by these equations 
is a signature of underdamped quantum phase dynamics. 
It occurs only if the impedance of the external circuitry 
is sufficiently large both at the frequency of Bloch oscil- 
lations and at the Josephson frequency of the individual 
elements. The precise conditions are given in Section fVl 
Observation of this dependence and the measurement of 
the maximal voltage would provide the proof of the quan- 
tum dynamics and the measurements of the tunneling 



amplitude which is the most important characteristics of 
these systems. It would also provide a crucial test of the 
quality of decoupling to the environment. 

As a deeply quantum mechanical system, the chain 
of Josephson devices is very sensitive to an additional 
ac driving. It exhibits resonances when the driving 
frequency is commensurate with the frequency lob = 
27r£7/2e of the Bloch oscillations. This would provide 
additional ways to characterize the quantum dynamics 
of these circuits and confirm the period doubling of the 
rhombi frustrated by exactly half flux quantum. 
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APPENDIX A: QUANTUM-MECHANICAL 
CALCULATION OF THE DC VOLTAGE 

In the large current regime where I > / max , the energy 
drop Ab = hi lie induced by the driving current when 
<f> increases by 2tt becomes comparable to or larger than 
the bandwidth W. In this regime, the semi-classical ap- 
proach is no longer reliable. But as long as Ab remains 
small compared to the typical gap A between nearby 
bands, we may still construct the system wave functions 
in the presence of the driving field from Wannier orbitals 
belonging to a single band. In such quantum-mechanical 
approach, dissipation is described as the result of cou- 
pling the single degree of freedom (0, q) to a continuum 
of oscillator modes (q a ,p a ). The corresponding Hamil- 
tonian has the form: 

H n =e n iq -Y<9 a q a I (Al) 

\ a / a 

where we have chosen the following commutation rela- 
tions: 



(A2) 



and all other commutators between these operators van- 
ish. The form of (|Aip is plausible on the physical ground 
because when the superconducting islands are coupled 
to macroscopic leads, the charge q undergoes quantum 
fluctuations, so that it has to be replaced by a "dressed 
charge" q — ^ Q g a q a in the effective Hamiltonian. A 
more explicit justification is that the corresponding semi- 
classical equations of motion have the same form as 
Eqs. (|4]), ([5]) which simply mean that the effective cur- 
rent going through the superconducting circuit is the bias 
current minus the current going through the external im- 
pedence. The semi-classical equations deduced from (|A1[) 
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read: 



# 


_ 1 de n 


dt 


h dq 


dq 


T 
1 


~dt 


Ye 


dq a 




dt 




dp a 




dt 


= -u a q, 



-(q -^2g a q a ) 



g a de n . ■sr-^ . 

It is then natural to introduce q' = q — ^ a g a Qot^ so that: 

(A3) 

E« 
3a- 



d$ _ 1 de n , 
dt h dq 
dq' I \r-^ dq a 



dt 2e 



dt 



(A4) 



To show that (|A4|) has the same form as (|3J|, we notice 
that the driving term for the bath oscillators is directly 
proportional to dtfi/dt. Specifically: 



dq a 2 d(p 
— +u; a q a =uj a9a - 



(A5) 



Going to Fourier space, we see that after averaging over 
initial conditions for the bath oscillators, Eq. (|A4[) takes 
the form: 

- iujq'(uj) = —2<k5(oj) - i ^ g 2 a ^1^2 

" a a 

(A6) 

This is exactly the frequency space version of Eq. |5]), 
where as usual, the dissipation is related to the spectral 
density of the bath by: 



uj 2 - uo 2 a Z(lu) 



(A7) 



Here we emphasize that Z is typically frequency depen- 
dent, in which case the term (Zq/ 'Z)(d<p/ 'dt) in Eq. ^ 
becomes a convolution with Zq/Z replaced by a non-local 
kernel in time. 

Now we turn to the solution of the quantum prob- 
lem (lAip in the large driving current regime. Let us first 
consider the Josephson array without dissipation. It is 
straightforward to express its eigenstates in the q repre- 
sentation^ because the Schrodinger equation reads then: 

E^( q ) = e n (q)TP(q)-i^(q) (A8) 

so that: 

1>(q) = 1>(0) exp (-i^ J\tnW) - E)dq^j (A9) 

The energy spectrum is determined via the boundary 
condition ip(q + 1) = ip(q) so that: 



.1/2 
1/ 



E v = e n (q')dq' + A ws i/, v integer (A10) 

.7-1/2 



This yields a Wannier-Stark ladder of spacially localized 
states, with a constant level spacing equal to Ab- Note 
that increasing v by one unit multiplies the wave- function 
ip(q) by exp(i27rg). In the phase representation, this is 
equivalent to a translation by —2tt. Of course, in the 
absence of dissipation, these levels have infinite life-time, 
and therefore, no dc voltage is generated. 

Let us now consider the limit of a weak coupling to 
the dissipative bath. This means that the decay rate T 
of the Wannier-Stark levels is much smaller than the level 
spacing Ab . Assuming that transitions take place mostly 
between two adjacent levels, we get an average voltage: 



(V) 



h ,d<j). 
2e 



hT 
~2e 



(All) 



The rate T is estimated via Fermi's golden rule which we 
prefer to use in the correlation function formulation : 

IW = h{v\^\v')\ 2 C AA {{v-v l )uj^) (A12) 
n A dq 

where lu-q = Ab/S = 2irl/(2e). In this expression, 
C aa is the Fourier transform of the correlation function, 
CAAit — t') = (A(t)A(t')) of the Heisenberg operators 
A(t) = ^2 Q g a q a (t), taken in the equilibrium state of the 
dissipative bath. 

We evaluate now the matrix element of the velocity 
operator de n /dq between Wannier-Stark states: 

1 /2 

H^V) = / ^r(«) exp(t27r(i/ - u)q) dq (A13) 
dq J-i/2 dq 

As we have seen, in most physically interesting situations, 
we can approximate the periodic function e(q) by a single 
harmonic 2wcos(27rg). In this case: 

(Hy>') = 2™V-H,i (A14) 

In the zero temperature limit, the bath correlation func- 
tion is: 

C AA {w) =J2^9 2 J{uj-uj a ) = - Re (^-) B{u) (A15) 



where 6(u>) is the Heaviside step function. Putting all 
these elements together gives: 



(V(I » J c )) = 



(2tt W ) 2 
2e 2 Z,J 



(A16) 



where Z u denotes the external impedance taken for 
w = bj-Q, and this result is in perfect agreement with the 
large current limit of the semi-classical treatment, shown 
as Eq. ([33]). 

When the band structure is replaced by 2w cos(27r£<7) 
as in the case of a rhombus (for which ( = 1/2), we have 
now: 



(v) = C 



hr 



(A17) 
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The frequency of Bloch oscillations becomes u>b(C) = 
C^b(C — 1)) an< i the matrix element is multiplied by 



C, so that the voltage is multiplied by C 2 - Again, this is 
compatible with the semi-classical result (I55|) . 
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